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1. Introduction 

Magnetic fields have been detected in planets, stars, galaxies, and clusters of galaxies. Intra¬ 
cluster magnetic fields arc characterized by a strength of ~ /rG and fluctuations on scales from a 
few kpc to ~ 100kpc (e.g., EnBlin & Vogt 2003, Vogt & EnBlin 2003, Vogt & EnBlin 2005, Murgia 
et al. 2004, Govoni et al. 2006, Guidetti et al. 2008, Bonafede et al. 2009, Bonafede et al. 2010, 
Vacca et al. 2010, Vacca et al. 2012). On larger scales, along galaxy filaments as well as in voids and 
sheets, a firm detection has not yet been possible. Cosmological magneto-hydro-dynamical sim¬ 
ulations indicate that magnetic fields are present in intergalactic low-density environments (e.g., 
Dolag et al. 2008), while y-ray observations suggest their presence in voids where a lower limit of 
~ 1 fG has been derived (see Beck A. et al. 2013 and references therein). For more details of our 
current observational understanding of magnetic fields in low-density regions, see also the chapter 
by Giovannini et al. (2015) in these proceedings. 

Magnetic fields can be studied using a variety of techniques, among which rotation measure 
(RM) observations are a powerful tool (e.g., Carilli & Taylor 2002, Govoni et al. 2004). The amount 
of Faraday rotation measured from radio observations along a given line of sight is the result of the 
contribution from the Milky Way, the radio source itself and other extragalactic environments. The 
distribution of matter in large-scale structures can be classified into four types (following Hahn 
et al. 2007, using a criterion of local stability of test particle orbits): voids, sheets, filaments and 
galaxy clusters. We will hereafter call these four types, as well as the sources responsible for the 
radio emission, environments. 

Disentangling the contributions from different extragalactic environments is essential in order 
to study magnetic fields in large-scale structures, and consequently infer information about cosmic 
magnetism (e.g.. Beck R. et al. 2013, Gaensler et al. 2004). However, this is a non-trivial task 
that demands a statistical approach. The first statistical studies to prove the existence of Faraday 
rotation from large-scale extragalactic environments (galaxy clusters) were published by Lawler & 
Dennison (1982) and Vallee et al. (1986). Later, an important statistical RM analysis on intracluster 
magnetic fields was performed by Clarke et al. (2001, 2004). Relevant contributions regarding 
intergalactic magnetic fields on even larger scales were done by Kolatt (1998), Stasyszyn et al. 
(2010), and Akahori et al. (2014). 

A statistical theory for Bayesian inference on spatially distributed signals. Information Field 
Theory, was developed by EnBlin et al. (2009). Within the framework of this theory, Oppermann 
et al. (2012, 2014) reconstructed the Galactic Faraday rotation foreground and estimated the ex¬ 
tragalactic contribution. In this chapter we take a step forward, with a new algorithm, able to 
statistically separate the different extragalactic contributions, by combining RMs with information 
about the large-scale structure of the cosmic web. This approach will make it possible to study 
magnetic fields on large scales with the RM grid data that will be obtained with the new generation 
of radio interferometers (see the chapter by Johnston-Hollitt et al. 2015 in these proceedings). 

2. General description 

The probability density distribution of the extragalactic Faraday rotation measure along a given 
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line of sight i can be described with a Gaussian of mean (0 e ,i) and variance (0 e 2 -) 


<0e,i> = 0 




J "e (*K (x') (B x (x)B x ( X ')) 
(1 + z(x)) 2 (l + z(x')) 2 


( 2 . 1 ) 


where a 0 is a constant, x, is the proper distance of the source, z is the redshift, n e is the electron 
density, and B x is the component of the magnetic field along the line of sight. We can define a length 
scale A*, so that the auto-correlation tensor of the magnetic field ( B x (x)B x (x')) on average gives 
significant contributions for \x — x'\ < A v , but not for \x — x'\ S> A v . The strength of the magnetic 
field, its correlation length, and the electron density have a different value in each environment j, 
i.e., in voids, sheets, filaments, galaxy clusters, and the in radio source itself. If we define 


Xx = aln t {x) 1 {B x {x) 2 )X x , 


( 2 . 2 ) 


which we assume to depend only on the environment, so that % x = Xj(x)> an d introduce fj as the 
length of the line of sight i through each environment j, the overall variance in Faraday rotation 
can be written as 



(x) 2 (B x (x) 2 ) A r 
(1 +z(x)) 4 


dx 


xo+YJuXj- 

j 


(2.3) 


The intrinsic contribution of the emitting source Xo is assumed to be on average the same for all 
sources. Since the intrinsic size of the emitting radio source is unknown, Xo is defined according to 
Eq. ( |2.2| ) but with the length /,o factored in. Here, the auto-correlation function of the magnetic field 
has been approximated with a top-hat function of width A v . We assume that, within a correlation 
length Ac, the redshift can be approximated to be constant. Moreover, we note that Xj has been 
assumed to be independent of the redshift, meaning that the contribution to the Faraday rotation 
from a given environment does not change with redshift. The redshift dependence is totally ab¬ 
sorbed by the term fj. Residual small-scale, uncorrelated Galactic contributions could be present 
in the estimate of the extragalactic Faraday rotation <j) c j. To properly take them into account, the 
inclusion of a latitude dependent contribution in Eq. (2.3) could be necessary. 

To constrain the contribution from the different environments Xj, an estimate of the extragalac¬ 
tic Faraday rotation (j) e j for a collection of lines of sight, as well as a description of the cosmic web 
large-scale structure, summarized by fj, is required. With a Bayesian approach, Jasche et al. (2010) 
reconstructed the cosmic density field on the basis of optical data from the SDSS Data Release 7. 
Following the scheme proposed by Hahn et al. (2007), they give a three-dimensional classification 
of the density field posterior in terms of voids, sheets, filaments, and clusters of galaxies. For a 
source with known redshift, this reconstruction allows one to estimate the path covered by the ra¬ 
dio signal through each environment, i.e., the elements fj in Eq. (2.3). Radio sources can sit in 
different structures, e.g. galaxy clusters or filaments. This effect can be taken into account statisti¬ 
cally since it is possible to locate the radio source in the reconstruction of the large scale structure 
via redshift identification. The fractions of the line of sight that the different radio sources spend in 
the different environments varies (this is also the case for the same radio source when different re¬ 
construction of the large scale structure are considered). The combination of this information with 
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the Faraday depth estimation for a collection of sources enables us therefore to statistically dis¬ 
criminate between the average contributions from the source itself, from galaxy clusters, filaments, 
voids, and sheets. 


3. Bayesian approach 


To infer information about the contribution from the different extragalactic environments, x> 
to the extragalactic Faraday rotation on the basis of data cl 1 obtained from a RM grid, we propose 
a Bayesian approach. Bayes theorem, 

P{d\x)Pix) 

P(d ) ’ 1 ; 

expresses our a posteriori knowledge P{%\d) about the signal % after the measurement process 
took place. The initial knowledge on the signal before the data have been acquired is described by 
the prior P(x) and is modified by the data d through the likelihood P(d\%). The evidence P(d) is 
a normalization factor, obtained by marginalizing the joint probability P(d,x) = P(d\x)P(x) over 
all the possible configurations of the signal X- 

In order to get a data-driven solution and to keep our assumptions as general as possible, 
we do not put constraints on the magnetic field strength and assume all values of magnetic field 
magnitudes to have equal probability. Such a priori knowledge can be well described by a flat prior 
on a logarithmic scale. Mathematically, it is expressed as an Inverse-Gamma distribution 

P(X)~ (3-2) 

j 

where cxj = 1 and c/j = 0 (Jeffreys’ Prior). 

The Faraday depths d, coming from radio observations contain the Galactic and extragalactic 
contribution, </)„., and </> e( -, and the noise /;, of the measurement process 

di — i T 0e,i H - tli. (3-3) 



Oppermann et al. (2012, 2014) developed an algorithm to reconstruct the Galactic Faraday rotation 
foreground and estimate the extragalactic contribution, given a catalog of Faraday RMs for different 
lines of sight. Their results are available as a collection of N random samples extracted from the 


posterior distribution 


P(<j>e|d,Opp) 


/V|<j>e)fQe|Opp) 

P(d |Opp) 


(3.4) 


which intrinsically depends on the assumption enclosed in their prior, as indicated by “Opp”. In 
order to disentangle the different contributions that make up the overall observed extragalactic 
Faraday depth, this prior must be replaced with a new one. 

’Here, d and % are respectively vectors with elements dj and %j with i = l,...,N\ os and j = 0, ...,N env , where N\ os 
and IVenv are the total number of lines of sight and the total number of environments in between the source and the 
observer. 
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If in Eq. ( p.l| ), we expand our likelihood 

P(d\x) = J ®&P(dmP(k\x) 


and insert the likelihood P(d\(j) e ) from Oppermann et al. (2014), see Eq. (|3~4|), we obtain 

P ( X= J ^P(^\d,Opp)P(d\Opp) P{&\x)P{x) 


P(0 e |Opp) 


P(d) 


(3.5) 


(3.6) 


where f &(j) c denotes the phase space integral over all the possible held configurations of 0 e . As 
the evidence is simply a normalization factor, this posterior no longer depends on the assumptions 
of Oppermann et al. (2014), which cancel out. 

Under the assumption that the lines of sight are sufficiently separated, that is, the cross¬ 
correlation function of then - magnetic held is zero, for each line of sight we can consider indepen¬ 
dent Gaussian distributions for the extragalactic Faraday depth given the cosmic web’s large-scale 
structure 

P{^\X) = I W(/|Jasche)n^(0 e ,i,£/ o ^), (3.7) 

where P(/|Jasche) is the posterior of the cosmic web structure from Jasche et al. (2010), available in 
the form of N' samples. The notation ( /j(x. : X) indicates a one-dimensional Gaussian distribution for 
a variable x with zero mean and variance X. Moreover, since the prior P(0 e |Opp) from Oppermann 
et al. (2014) is 

/ , (0 e |Opp) = n^eAtfe) (3-8) 

i 

with a standard deviation <7 e ~6.4 rad/m 2 , it follows that 


n y^ii% iij Xi) 

P(X,$tAd) ~ ? ricefA 1 2 \ n^ 7lp (^l^’ Q PP) /3 ( / l Jasche )- ( 3 -9) 

I | (0e,(, <J C ) j 

i 


This probability density function can be marginalized over 1 and 0 e by making use of the sample- 
representation available for the posterior of the large scale structure P(71 Jasche) and of the extra- 
galactic Faraday rotation P(0 e |d,Opp). To evaluate the posterior distribution for we plan to use 
a Monte Carlo Markov Chain approach. 

In this approximation the lines of sight are considered independent. This means that, if the 
assumption does not hold, a correlated component of the magnetic held could be missed. A proper 
modeling that takes into account possible cross-correlations between lines of sight is left for future 
works. 


4. Current developments and applications with SKA 

An application to archival Very Large Array (VLA) data is in progress. By reanalyzing 1.4 GHz 
data from the NRAO VLA Sky Survey (NVSS), Taylor et al. (2009) derived a catalog of RM values 
for 37543 linearly polarized radio sources, north of declination —40° - which corresponds to a 
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number density of approximately one source per square degree. Oppermann et al. (2014) evaluated 
the posterior distribution of the extragalactic Faraday rotation for all the lines of sight in this catalog. 

By using optical data from on-line databases (NED and SIMBAD) as well as from surveys 
(SDSS, 6dFGS, 2dFGRS, 2QZ, 6QZ), Hammond et al. (2012) published a catalog of spectroscopic 
redshifts for 4003 linearly polarized sources from the NVSS. The sources in this sample have 
redshifts in the range 0.0007 < z < 5.3. Only about 30 sources have z < 0.02, and the mean 
redshift of the sample is 0.89. The mean angular distance of the sources being equal to ~1 degree, 
this implies a mean distance among sources larger than 1 Mpc for 99% of the sample, and on 
average larger than ^30 Mpc. If A v >>30 Mpc, the assumption of no cross-correlation is not valid 
anymore. As noted in the previous section, this would translate in an underestimate of the magnetic 
held for structures with such a length scale. By combining the reconstruction of the large-scale 
structure from Jasche et al. (2010) with this catalog of redshifts, it is possible to evaluate the path 
covered by the radio signal through each environment for each line of sight. 

The catalog of Faraday rotation values in Taylor et al. (2009) was derived through a linear ht 
of the observed polarization angle versus A 2 at two closely separated radio frequencies (1365 and 
1435 MHz). When few data points are used, such A 2 -hts can suffer from an n ^-ambiguity. More¬ 
over, this approach is suitable only if no Faraday rotating plasma is mixed with the radio-emitting 
plasma. A better analysis of polarization properties of radio sources can be done with the RM 
synthesis (Brentjens & de Bruyn 2005) and the Faraday synthesis (Bell & EnBlin 2012) techniques 
if more frequencies are available. A continuous and broad frequency coverage is desirable, both 
to reduce the risk of //Tt-ambiguitics and to reach a sufficiently high resolution in Faraday depth to 
distinguish nearby Faraday components (e.g., Farnsworth et al. 2011, Kumazaki et al. 2014). 

In this respect, the Square Kilometer Array as well as its precursors will offer a great improve¬ 
ment. Broadband spectro-polarimetric surveys of the sky in combination with deep observations 
of specific targets will provide densely spaced RM grids. A RM catalog one hundred times denser 
than those currently existing (Taylor et al. 2009) will be obtained with POSSUM, the polarization 
survey planned with the SKA precursor ASKAP. POSSUM will observe the sky between 1130 
and 1430 MHz with a sensitivity in U and Q of ~ 10 /r Jy/bcam and a resolution of 10", giving an 
expected density of 70 polarized sources per square degree (Hales et al. 2014a). The uncertainty 
on the RM is lower than ~ 6 rad/m 2 for a polarized signal with a S/N >5 (Stepanov et al. 2008). 
The ultra-deep polarization survey CHIFES Con Pol with the Jansky Very Large Array (Hales et al. 
2014b) will reach a sensitivity of 400 nJy with sky area of a fraction of a square degree, providing a 
first look at the magnetic held science that can be performed with the SKA. The sky survey in polar¬ 
ization SKA 1-SUR planned for SKA1 will go much further. Designed to cover the frequency range 
650-1670 MHz, it will produce Stokes U and Q images over 30000 square degrees with a spatial 
resolution of 2". Within two years, a sensitivity in polarization of ~ I /./Jy/bcam will be reached, 
with an uncertainty in RM lower than ~ 1 rad/m 2 (.S’ /N > 5). The high RM accuracy will make our 
approach viable for detecting magnetic fields with strengths ~ nG already during the early science 
phase of SKA1, when sensitivity will be at 50% of its full level. For weaker magnetic fields the full 
deployment of all the SKA1 baseline will be necessary. The resulting sky grid of Faraday rotation 
values will be 300-1000 times denser than the largest catalog available at the moment. The better 
resolution of 2", compared to the 45" of Taylor et al. (2009), will make it possible to identify optical 
counterparts uniquely and hence the source redshift. For fields of particular interest, e.g., regions 
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where galaxy filaments have been observed, deeper observations could be necessary to assess the 
presence of magnetic fields and their properties. In this case, the observations from the survey can 
be complemented by targeted SKA 1-MID pointings over a continuous frequency band from 350 to 
3050 MHz. This would imply a polarization sensitivity < 1 /./Jy/bcam after one hour of integration 
time and an uncertainty in RM lower than ps 0.2 rad/m 2 ( S/N > 5). These targeted observations 
will be crucial, both to increase the resolution in Faraday depth and to detect a larger number of 
sources, allowing to distinguish nearby components in Faraday spectra and improve the statistics. 

The even higher spatial and frequency resolution and the deeper sensitivity of SKA2 will make 
it possible to obtain RMs for an unprecedented number of background and intervening sources. It 
is therefore essential to develop and apply statistical approaches with the aim of interpreting the 
data correctly. 

5. Conclusions 

The origin and evolution of cosmological magnetic fields arc still open questions. Magnetic 
fields arc currently being studied in galaxy clusters but a firm detection is lacking in larger scale 
structures, such as galaxy filaments. The structure and the strength of magnetic fields in galaxy 
clusters arc both deeply affected by the cluster formation processes (e.g., Dolag et al. 2005, Don- 
nert et al. 2009, Xu et al. 2010), almost independently of the mechanism of generation of the seed 
magnetic field. This makes it difficult to discriminate between the different formation scenarios 
from galaxy cluster observations alone. The SKA’s large bandwidth, high sensitivity and frequency 
resolution will permit a more reliable recovery of the Faraday depth for a large number of sources, 
and its high spatial resolution will allow a proper identification of the optical counterparts neces¬ 
sary for redshift assignment. Modern techniques based on the largest photometric all-sky surveys 
(2MASS, WISE and SuperCOSMOS) already provide catalogs of redshifts for millions of sources 
up to redshifts of 0.2 to 0.3 (e.g. Bilicki et al. 2014), and the next generation of redshift surveys 
will push those limits. We will investigate to what depths and accuracy redshifts are needed in 
order to apply the proposed method. The ultimate goal is to search for magnetic fields in large- 
scale structures like filaments, where magnetic field amplification is not yet saturated and therefore 
the field strength still depends on the seed field strength. This information will help to place firm 
constraints on the origin and evolution of cosmological magnetic fields. 
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